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NOTATION 

Dotted  variables  are  dimensional j  undotted  variables  are 
dimensionless 


a,  b 


f,  P 


PrT  =  U,2/gT ' 
PrL  =  U'2/gL' 


constants 
ship  beam 
drag  force 

outer  and  inner  complex  potentials 


Froude  numbers 


g  gravity  acceleration 

h'(x',z'),  h'(x')  functions  describing  the  hull  shape 


k(z) 


U,  V,  W 


function  of  complex  variable  (k  =  w  +  if) 

shiplength 

forebody  length 

free  surface  elevation  (inner,  dimensionless) 
pressure 

pressure  (inner,  dimensionless) 

velocity  modulus 

draft 

velocity  components 

velocity  component  (Sect.  II),  complex 
velocity  w’  «  u'  -  lv'  (Sect.  Ill) 

unperturbed  velocity  at  infinity 

velocity  components  (inner  dimensionless 
Sect.  Il)j  w  ■  u  -  iv  (Sect.  Ill) 

velocity  vector 

coordinates 
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z ' 

X,  Y,  Z 


coordinate  (Sect.  II),  complex  variable 
z1  =  x'  +  iy'  (Sect.  Ill) 

coordinates  (inner,  dimensionless) 


y 

A' 

€  =  T/L‘ 

eB  =  B,/L' 

€*  =  T'/gU'2 

<t> 

(D,  C 
i,  M- 
X 
9 

n  —in  (U-iV) 


constant 

jet  thickness  (A  =  A'/T') 
draft/length  ratio 
beam/length  ratio 

velocity  potential 

velocity  potential  (inner  dimensionless) 

auxiliary  complex  variables 

auxiliary  variables 

angle  and  also  dummy  variable 

angle  between  velocity  vector  and  x  axis 
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I .  INTRODUCTION 

The  wave  pattern  created  by  a  ship  moving  steadily  in  an 
ideal  fluid  and  the  related  wave  resistance  are  classical  subjects 
of  hydrodynamics.  Although  the  theory  has  diversified  and  com¬ 
putational  refinements  have  been  achieved  with  time,  there  has 
been  little  essential  progress  beyond  the  linearized  techniques 
introduced  by  Michell  and  Havelock.  In  essence  the  present  method 
of  solution  of  the  gravity  flow  problem  is  based  on  two  approxi¬ 
mations:  (i)  the  free  surface  condition  is  linearized  and  (ii) 

the  hull  is  replaced  by  a  singularity  distribution  along  a  line 
or  a  plane.  The  wave  resistance  is  generally  determined  from  the 
rate  of  energy  radiated  far  away  from  the  ship.. 

The  above  two  basic  approximations  have  been  given  foundation 
in  a  rational  way  in  the  last  years  by  the  application  of  the 
method  of  matched  asymptotic  expansions  (Tuck  1965*  Ogilvie  1967). 
It  has  been  shown  that  the  classical  theory  is  in  fact  a  first 
order  term  of  an  outer  expansion  in  which  the  observer  is  fixed 
with  respect  to  the  ship  length  while  the  ratio  draft/length  or 
beam/length  (or  both)  tend  to  zero  and  the  Fr  number  based  on 
length  remains  constant.  In  the  vicinity  of  the  body,  in  the 
inner  zone,  the  solution  is  still  valid,  provided  that  the  slen¬ 
derness  parameter  is  sufficiently  small  and  the  ship  has  a  fine 
form.  The  latter  condition  has  been  somehow  overlooked  when 
applying  the  theory  to  actual  ships  which  do  not  generally  have 
a  needle-like  or  knife-like  shape.  In  the  extreme  case  of  a 
blunt  shape  there  is  stagnation  at  the  bow  and  the  linearized 
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assumptions  are  badly  violated  there:  the  speed  and  the  free 
surface  rise  are  no  longer  small  perturbations  of  the  uniform 
speed  and  of  the  horizontal  level,  respectively.  At  the  stern 
the  situation  is  different  due  to  separation  and  viscous  effects. 

Realizing  the  importance  of  bluntness  effects  on  resistance 
of  real  ships,  we  have  initiated  a  study  of  the  free  surface  flow 
near  the  bow  and  of  the  related  resistance. 

The  present  report  summarizes  our  first  results  which,  be¬ 
cause  of  the  complexity  of  the  problem,  involve  in  thi3  initial 
stage  rather  crude  approximations  obtained  for  highly  schematized 
configurations.  It  is  our  feeling,  however,  that  this  initial 
step  yields  basic  understanding  of  the  problem.  We  hope  to  be 
able  to  extend  and  refine  the  results  in  the  future,  to  compare 
them  with  experiments,  and  eventually  to  apply  them  to  actual 
ships . 

II.  INNER  AND  OUTER  EXPANSIONS  AND  CLASSIFICATION  OP  SHIPS 
1.  Notation  and  Basic  Equations 

The  symbols  used  in  this  report  arc  given  in  "Notations" 
and  also  shown  in  part  in  Figure  1. 

The  ship  le  defined  by  the  shape  of  its  hull  represented 
by  the  equation 

f(x',  y*,  *' )  -  0 
or  in  one  of  the  explicit  forms 


[2.1] 
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y*  =  h'(x',  z')  [2.2] 

The  three  basic  lengths  associated  with  the  hull  are  B',  T* 
and  L'.  Additional  geometrical  coefficients  or  lengths  may  be 
considered,  like  the  forebody  length  Lx'.  The  latter  is  im¬ 
portant  in  characterizing  the  bluntness. 

Assuming  that  the  flow  is  steady  and  uniform  at  infinity, 
the  equations  satisifed  by  V'  and  tj  • ,  given  here  for  convenience 
of  reference,  are  as  follows: 


rot  V'  «  0  * 

(2.3 

div  V'  -  0 

J  (in  the  flow  domain) 

[2.4] 

uf*+  v,a+  w,a 

U'a 

(2.5) 

2  +  8ri  ■ 

2 

(on  the  free-surface 

>  y'  ■  y  (x * , z 1 ) ) 

v-  u'y,x,  -  wy,z, 

-  C 

J 

I 

(2.6] 

v  -  u'h* 

'x«  ‘  *  0 

1 2-71 

Equations  [2.3)  and  (2.4]  express  as  usual  irro tat Iona llty 
and  incompressibility,  Equation  (2.5)  is  the  dynamic  Bernoulli 
condition  on  the  free  surface,  while  Equations  [2.6]  and  (2.7) 
are  the  klnematlcel  boundary  conditions  along  the  free-surface 
and  the  hull  respectively. 
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In  order  to  render  the  solution  unique,  the  radiation  con¬ 
dition  is  imposed 


U1  =  -U';v'  =  W'=0  (x1  -*■  +  oo  j  y<  -*•  -  oo )  [2.8] 

Equations  [2.3]  to  [2.8]  may  be  reformulated  in  terms  of  the 
velocity  potential  <t> 1  by  replacing  V  by  grad  <!>'. 

2.  Outer  Expansion;  Classification  of  Hull  Shapes 

Economy  is  achieved  by  making  the  variables  and  the  equa¬ 
tions  dimensionless  in  the  standard  way.  Let  us  define  the 
following  outer  variables 

V  -  V'/U'(u,v,w  =  u’/tJ',v'/U',w'/tJ' );  x,y,  z  =  x '/L' ,  y '/L' ,  z ; 

r)  =  V/L'J  h  ■  h'/L';  p  =  p'/pU,2;<t>  =  <D'/U'L' 

[2.91 

The  equations  of  flow  [2.3]  -  [2.8]  become,  in  terms  of  <t>. 


Vz<fc  -  0 


(in  the  flow  domain) 


r\/Frrz  «  1/2 


>  (y  -  n(x,z)) 


’.y  ■  *'xn'x  ■  *’2^2  ’  °J 


*>y  -  -  *>zh>z  "  0  (y  *  h(x,z)) 


-  -  *  ) 


[2.10] 

[2.11] 

[2.12] 

[2.13] 


u  =  -1  j  v  ■  w  ■  0 


(x  -*•  +  00  j  y 


[2.14] 
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T'he  solution  of  <t>  depends  on  x,  y,  z  and  the  parameters 

x 

PrL  =  )2>  B'/L'j  T'/L'  for  a  given  hull  shape.  The  clas¬ 

sical  technique  for  simplifying  the  nonlinear  problem  is  to  take 
advantage  of  the  fact  that  B'/L'  or  T'/L'  or  both  are  much  smaller 
than  unitjr.  With  e  =  T'/L'  and  =  B'/L'  an  outer  expansion  is 
obtained  by  assuming  that  V  and  $  may  be  expressed  as  a  series 

associated  with  an  asymptotic  sequence  based  on  e  or  €„.  This 

B 

has  been  done  in  numerous  publications  (see  for  instance  Wehausen 
and  Laitone,  i960  and  Tuck,  1965)  and  will  not  be  repeated  here. 
Since  by  definition 


h  (x,z)  =  e(x,z)  [2.15] 

where  H  =  0(l),  it  is  natural  to  consider  an  expansion  of  the 
type 

<t>  =  -  x  +  ed>i  +  0(e ) 

T]  =  €T|1  +  0(e)  [2.16] 

with  x,  y,  z  ■  0(l).  For  an  outer  observer  with  a  position  fixed 
with  respect  to  the  ship  length,  the  ship  collapses  in  a  line  or 
a  plane  at  zero  order  and  the  flow  Is  unperturbed.  At  first 
order  (and  we  consider  here  only  first  order  terms)  the  equations 
become  the  well  known  linearized  equations 
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V2<J>1  =  0 

y  <  0 

[2.17] 

-ui  +  T)i/FrL2  =  0  *j 

j 

>  (y  -  0) 

[2.18] 

Vl  -  T)i,x  =  0  J 

[2.19] 

Vl  +  H,x  =  0 

(on  the  hull) 

[2.20] 

u  =  vi  =  Wi  =  0 

(x  -*  +  00  ;y  -*•  -  °°) 

[2.21] 

The  hull,  at  first  order,  degenerates  into: 

(  i)  a  line  in  the  case  of  slender  ships  (eB  =  0(e)), 
(ii)  a  vertical  plane  in  the  case  of  thin  ships 
(e  =  0(l )  and  e0  =  0 ( 1 ) )  and 
(iii)  a  horizontal  plane  at  y  =  0  for  flat  ships 
(efi  =  0(1)  and  e  =  0(l)). 

Different  flow  regimes,  and  equations  accordingly,  are  ob¬ 
tained  corresponding  to  the  relationship  between  Pr^  and  e. 
Ogllvie  (1967)  has  analyzed  these  possibilities.  Since  we  con¬ 
centrate  here  on  displacement  ships  mainly,  we  should  consider 
the  following  two  possibilities: 


(i)  Small  PrL  number,  Pr^2  =  0(l).  In  this  case  a 

direct  expansion  of  Equations  [2.10]  -  [2.11]  gives 


vi  =  0 


(y  =  0) 


[2.22] 
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i.e.,  a  rigid  wall  condition  of  the  free  surface.  This  expansion 
is  discussed  in  detail  in  Section  III. 3. 

(ii)  FrL  =  0(l),  which  yields  the  ordinary  ship  re¬ 
sistance  problem,  with  gravity  waves  left  behind 
the  ship  (Equations  [2.17]  -  [2.21]) 

Higher  Fr^  lead  to  planing  problems  not  considered  here. 

Problem  (ii),  by  far  the  most  interesting,  has  been  solved 
by  replacing  the  degenerated  hull  by:  (i)  a  line  of  sources  for 
slender  ships,  (ii)  a  source  distribution  in  the  mid-plane  for 
thin  ships  and  (iii)  a  pressure  distribution  on  the  free  surface 
for  flat  ships  (Lunde,  1952). 

The  fulfillment  of  the  free-surface  conditions  [2.18]  and 
[2.19]  is  equivalent  to  the  extension  of  the  flow  in  the  whole 
space  (above  y  =  0)  and  the  introduction  of  an  infinite  system  of 
singularities  in  y  >  0,  reflection  of  the  ship  singularities.  In 
the  case  of  small  FrL  just  one  image  is  sufficient  in  order  to 
satisfy  Equation  [2.22], 

3 .  Inner  Expansion:  Bow  Singularity 

The  outer  expansion  is  singular  near  the  body  in  the 
FrL  =  0(l)  case  since  the  first  order  velocity  tends  there  to 
infinity.  Tuck  (1965)  has  considered  an  inner  expansion  for 
slender  ships.  The  inner  variables  are  defined  as 

X  =  x  -  x'/L';  Y  =  y/e  =  y'/T'J  Z  =  z/e  *  z'/T'J 


N  =  T]/€  =  rj'A'j  Uj  V,  W  -  u,  v,  w 


[2.23] 
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The  equations  of  flow  [2.10]  and  [2.11]  are  again  expanded 
by  assuming  that  U,V,W  and  N  are  asymptotic  series  with  respect 
to  e 


U,V,W  =  U  ,V  ,W  +  0(e) 
0  0  0 

N  =  N  +  0(e) 
o 


[2.24] 


while  X,Y,Z  =  0(l). 

In  the  inner  limit  the  observer  is  fixed  with  respect  to  the 
beam  (or  draft).  For  such  an  observer,  when  e  -*■  0  the  hull  cross 
section  keeps  its  shape  unchanged  while  the  shiplength  tends  to 
infinity.  The  equations  become,  at  zero  order,  two-dimensional 
(in  Y, Z)  and  the  free  surface  condition  becomes  that  of  a  rigid 
wall. 

The  matching  of  the  outer  and  inner  expansions  yields  (Tuck, 
1965)  the  classical  result,  at  first  order  :  i.e.  the  replacement 
of  the  ship  by  a  source  system.  The  source  strength  is  propor¬ 
tional  to  the  cross-section  area  variation.  The  inner  expansion 
is  valid  only  if  this  variation  is  gradual,  i.e.  for  fine  ships. 
The  slender  body  expansion  fails  in  the  bow  region  if  the  ship 
has  some  bluntness,  and  there  the  inner  problem  is  no  more  one 
of  two-dimensional  flow  in  the  Y, Z  plane,  nor  is  the  condition 
on  the  free  surface  one  of  a  rigid  wall.  For  this  reason  we 
should  call  the  slender  body  expansions  outer  and  inner  midbody 
expansions,  in  order  to  stress  their  limitations. 
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In  the  case  of  a  blunt-bow  ship  the  appropriate  inner  vari¬ 
ables  in  the  bow  region  are 

X,Y,Z  =  x/e,  y/e,  z/e;  U,V,W  -  u.VjW;  N  =  r\/e;  $  »  <J>/e 


U,V,W  =  Uo,Vo,Wo  +  0(e);  *  =  4>Q(X,Y,Z)  +  0(e)jN  =  N  +  0(e) 


[2.25] 


[2.26] 


with  X, Y, Z  =  0(1). 

While  in  the  case  of  the  inner  midbody  expansion  the  ob¬ 
server  is  fixed  laterally  with  respect  to  the  ship  and  at  zero 
order  the  length  tends  to  infinity,  in  both  bow  and  stern  direc¬ 
tions,  in  the  bow  inner  expansion  the  observer  is  fixed  with  re¬ 
spect  to  the  bow  and  the  shiplength  tends  to  infinity  sternwise. 

Substituting  [2.25]  and  [2.26]  into  Equations  [2.10]  - 
[2.13]  we  get  at  zero  order 


Va<t>Q  =  0  (in  the  flow  domain) 


[2.27] 


Ua+Va+Wa  ,  . 

0  0  0  1  „  1 

111  ■■■--'  1 - +  -  N  ■  — 

2  a  0  2 


V„  -  U  N  _  -  W  K  .  -  0 

0  0  0,X  o  0,2 


(Z-No) 


[2.28] 


[2.29] 


Vo  -  UoH’x  -  WoH'z  -  0 


(Z  -  H)  [2.30] 
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The  condition  at  infinity  [2.14]  is  lost  and  is  replaced  by 
the  matching  with  the  outer  expansion. 

i 

The  bow  Proude  number  Frm  =  U'/(gT')2  is  related  to  Fr_  and 

1  1j 

e  through 

FrT2  =  Fr/Ve  [2.31] 

Consequently  the  Bernoulli  Equation  [2.28]  may  have  the  following 

form  depending  on  the  order  of  magnitude  of  Fr_ : 

Li 

(i)  FrL2  =  0(l)(l/FrT2  =  0(e)).  Then  Equation  [2.28] 

U2+V2+W2»l  (Z  =  N  )  [2.32] 

ooo  'o' 

i.e.  free  gravity  flow  at  zero  order  in  the  inner  region. 

(li)  FrL2  ■  0(e)(l/FrT2  =  0(l)).  Equation  [2.28]  remains 

unchanged  and  we  have  the  full  nonlinear  gravity  problem, 
(iii)  Fr^2  =  0(e2 )(l/FrTa  •  0(l/e)).  This  case  reduces  to 
that  of  a  rigid  wall  condition  discussed  in  the  pre¬ 
ceding  section. 

In  the  case  of  slender  and  thin  ships  the  flow  in  the 
vicinity  of  the  bow  is  three-dimensional  in  all  the  above  approxi¬ 
mations.  Further  simplifications  are  achieved  in  the  case  of 
flat  ships.  Then  the  proper  inner  variables  are 


X  -  x/s;  Y  -  y/cj  Z  ■  z;  N  »  tj/c;  $  ■  ♦/*;  U*V,W  ■  u,v,w  (2.33) 
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while  the  asymptotic  expansion  starts  as 

<t  =  <Do  +  0(e);  N  =  No  +  0(e);  U,V  =  UQ,Vo  +  0(e);  W  -  0(e) 

[2.3*1 

Again  the  substitution  of  [2.33]  and  [2.3*1  Into  Equations  [2.10] 
-  [2.14]  gives,  at  zero  order 

[2.35) 

[2.36) 

(Y  -  No) 

[2.37) 

(Y  -  H)  [2.38] 

and  the  problem  is  reduced  to  that  of  gravity  flow  in  the  ver¬ 
tical  X, Y  plane  in  the  vicinity  of  a  body  of  shape  Y  ■  H(X,Z) 
(here  Z  appears  as  a  parameter).  The  requirement  of  flatness 
has  the  meaning  of  T'/B'  <  1,  but  still  allows  for  B'/L'  <  1. 

In  the  flat  ship  approximation  the  observer  attached  to  the  bow 
sees  both  width  and  length  tending  to  infinity  (although  pos¬ 
sibly  at  different  rates).  Obviously  this  approximation  is  not 
valid  near  corners  or  regions  of  large  change  of  H  with  X.  There 
the  full  three-dimensional  flow  or  some  other  approximations  have 
to  be  considered.  Again  we  obtain  the  three  different  cases 
discussed  above  depending  on  whether  l/Pr^**  0(c), 0(1)  or  0(l/t). 


7  <*>  «  0 

x,y  o 


°o“+V  +  1  H  i 
— r- 


V  -  U  N  -  0 
o  0  0, X 


v0  -  u0Hlj(  -  0 
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We  summarize  the  discussion  of  all  the  encountered  cases 
in  the  following  table: 


DISPLACEMENT  SHIPS  (PrT  <  l) 

J-l 


0(1) 

0(1) 

PrL  =U '  /gL 

0(e2) 

0(e) 

wBm 

0(e) 

0(1) 

0(e-1) 

Outer 

Expansion 

A 

Rigid  Wall 
Condition 

Linearized  gravity  waves  far 

« 

Rigj 

Cone 

.d  Wall 
lition 

from  the  ship.  The  ship  is 
replaced  by  a  line  (slender) 
or  a  plane  (thin  and  flat) 
distribution  of  singularities. 

Inner 

Midbody 

Expansion 

Everywhere 
The  un¬ 
perturbed 
flow  is 
the  state 
of  rest. 

Rigid  Wall 
Condition 

Rigid  Wall  Condition 

Slender  ships:  two-dim.  flow 
in  vertical  planes  normal  to 
the  centerline  (Tuck,  1965). 

Inner 

Bow 

Expansion 

mummsam 

Nonlinear  free-gravity  flow 

Flow.  Slender  and 
thin  ship?:  three- 
dim.  flew;  flat 
ships:  two -dim. 
flow  in  vertical 
planes  normal  to 
the  bow. 

Slender  ships:  three-dim. flow. 

Thin  ships:  three-dim.  flow 
near  a  strut. 

Flat  ship:  two-dim.  flow  in 
vertical  planes 
normal  to  a  body 
of  infinite  length. 

obs.:  In  the  case  of  fine  ships  the  midbody  expansion  is  valid 
everywhere. 
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Finally  we  present  in  Figure  2  a  plot  of  B'/T'  and  FrT  for 
wore  than  one  hundred  existing  ships. 

There  is  no  apparent  correlation  between  the  two  parameters. 
At  any  rate  most  of  the  ships  considered  are  flat  rather  than 
thin  (B'/T'  *  2.2  +  3.4). 

The  draft  Froude  number  Fr^  is  of  order  one  in  most  cases, 
but  reaches  values  as  high  as  2  for  a  rapid  containership  and 
more  than  3  for  cruisers  and  destroyers. 

III.  GRAVITY  FLOW  PAST  TWO-DIMENSIONAL  BLUNT  BODIES  OF 

SEMI-INFINITE  LENGTH 

1.  General 

In  the  preceding  sections  it  was  shown  that  in  the  case  of 
flat  ships  the  inner  bow  flow  reduces  to  a  two-dimensional  flow 
in  a  vertical  plane  normal  to  the  bow.  In  the  remaining  sec¬ 
tions  of  this  work  we  consider  exclusively  such  flows.  More¬ 
over,  we  are  assuming  that  the  outer  flow  is  also  two-dimensional 
and  that  the  body  is  of  semi-infinite  length.  Obviously,  these 
assumptions  simplify  the  problem  considerably.  The  essential 
features  of  the  bow  flow  are,  nevertheless.  Included  in  the 
picture.  We  plan  to  apply  the  results  by  some  approximate  tech¬ 
niques  to  actual  ships  in  the  future,  taking  advantage  of  the 
fact  that  for  most  ships  the  ratio  draft/beam  is  smaller  than 
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Conslstent  with  the  range  of  PrT  considered,  which  apply 
to  displacement  ships,  we  assume  that  the  bottom  of  the  midbody 
is  horizontal.  The  results  permit  however,  to  compute  trim  to  a 
first  approximation,  but  we  do  not  consider  this  problem  here. 

In  the  case  of  a  two-dimensional  flow  (Figure  3)  the  dimen¬ 
sionless  velocity  potential  depends  on  only  one  parameter  for 
a  given  hull  shape:  ♦  *  *  (x,y;FrT).  Consequently,  the  possible 
asymptotic  expansions  of  the  exact  Equations  [2.10]  -  [2.14], 
with  the  z  components  deleted,  reduce  to  the  following  cases: 

(i)  Small  FrT  .  In  this  case  the  Bernoulli  equation 
gives  a  rigid  wall  condition  in  a  first  approximation.  A  uni¬ 
form  expansion  solves  the  problem.  Results  for  the  first  and 
second  order  approximations  are  given  in  Section  III. 3. 

(ii)  Large  Fr^.  In  this  case  the  outer  flow  conforms 
to  equations  similar  to  the  linearized  Equations  [2.17] -[2. 21], 
while  the  inner  flow  is  that  of  a  free-gravlty  flow  at  zero 
order.  For  this  regime  we  sugg  s t  two  possible  inner  flow 
models:  The  Jet  model  discussed  in  detail  in  Section  III. 4  and 
the  spiral  vortex  model.  It  is  presumed  that  the  Jet  model  is 
adequate  for  large  Pr^,  while  the  spiral  vortex  model  represents 
moderate  to  large  PrT  flows.  Only  the  theory  for  the  former  is 
presented  herein. 

In  addition  a  discussion  of  the  exact  equations  of  free- 
surface  gravity  flow  near  a  stagnation  point  is  presented  in 
Section  III. 2. 
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2.  Free-Surface  Gravity  Flow  near  a  Stagnation  Point 

Let  us  consider  the  confluence  between  a  free-surface  and  a 
rigid  wall  in  the  vicinity  of  a  stagnation  point  (Figure  4a). 

In  the  symmetrical  case  (Xi  *  tt  -  12)  the  classical  Stokes 
result  (Wehausen  and  Laitone,  i960)  requires  that  X  *  X2  -  Xi 
=  120°.  This  result  will  now  be  extended  for  other  possible 
angles  between  AO  and  OB. 

In  the  vicinity  of  0(z  =  0,  Figure  5a)  we  assume  that  the 
z-plane  is  mapped  on  the  complex  potential  pla^  f  (Figure  5b) 
by 

-iXx  X/ir 

z  «  a  e  f  +  R(f )  [3.1] 

AOB  being  obviously  a  stream  line. 

The  function  R,  which  has  to  vanish  at  0,  Is  assumed  to  be 
in  the  vicinity  of  0  of  the  form 

R  -  b  fY  [3.2) 

wltn  b  -  b'ei0  a  complex  number  and  y  a  real  number. 

Obviously,  y  >  X/r,  otherwise  the  mapping  of  the  comer 
AOB  is  not  ensured. 

In  order  to  apply  the  Bernoulli  equation  along  AO  let  us 
determine  y  and  q*  -  u*  ♦  v*  as  functions  of  ♦.  From  Equa¬ 
tions  t3*l]  end  13.2]  we  obtain  on  AO  (f  *  e) 

X/r 

+  b'  sin  6  e 


y  *  -  a  sin  X*  ♦ 


(3.31 
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l/Q2  -  x,/  +  y =  (2lj\.<V*-i)  +  2  ab^x  oos(^I+6)» 


By  expanding  Equation  .[3-4],  Q2  is  found  as 


{\/r-y-2) 


[3.4] 


r  \2  -2(x/7r-l)r  Trb 1  y  Y -X/ir  1  r 

q2  _  (j,  1-2  — r-*-  0  cos(Xi+6)  +. . .  .  [3.5] 

Xa  a  A 


Substituting  y  and  q2  into  Bernoulli's  equation  and  retaining  terns 

r*  4*  *2  ~  4-  mnof  ua  rria4" 


of  order  0 


or  0  at  most  we  get 


a2  X/tt  7  1  /  ir\2  -2(X/tt-1) 

y  +  2g  “  "  3  sln  ^1<t>  +  13 '  Sin  54,r  +  2gXa’J<,>  +...=  0 

[3.6] 

The  identity  [3.6]  yields  the  following  relationships  be¬ 
tween  Xi  and  X: 

(i)  Xx  ^  0.  The  first  two  terms  of  [3.6]  give 


X  =  2tt/3 


§  —  =  g  sin  Xi 
8  a3 


[3.7] 


This  is  Stokes  classical  result.  Obviously,  X2  >  2tt/3 

(li)  Xi  =  0.  The  first  term  of  Equation  [3.6] 
vanishes,  and  the  remaining  give 


Y  =  -  2  (  X/tt  -  l) 


[3.8] 
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Slnce  y  >  X/tt  Equation  [3-8]  shows  that  X  <  2tt/3.  A  particular 
case  is  that  of  y  =  1,  which  renders  the  function  R  analytical. 

In  this  case  Equation  [3-8]  gives  X  =  7r/2,  i.e.  the  confluence 
between  a  horizontal  free-surface  and  a  vertical  wall. 

In  conclusion  there  are  two  possible  angles  between  a  free- 
surface  and  a  rigid  wall  at  the  stagnation  point: 

(i)  if  the  wall  is  inclined  with  respect  to  the  hori¬ 
zontal  at  an  angle  larger  than  120°(2tt/3  <  X2  <  t r) 
the  free  surface  intersects  the  wall  at  120° 

(X  =  2tt/3)  and, 

(ii)  if  the  wall  is  inclined  at  less  than  120° 

(X  <  2ir/3)  the  free-surface  is  horizontal 
(Xx  =  0). 

We  will  consider  blunt  bows  of  the  latter  type  in  Sec¬ 
tion  III. 3. 

3.  Small  Froude  Number  Flow  (Fr^,  <  l) 

(a )  General 

We  consider  here  the  flow  past  a  blunt  body  of  the 
shape  of  Figure  5a.  An  asymptotic  expansion  with  Fr,p  as  a  small 
parameter  has  as  its  zero  order  term  the  state  of  rest  cor¬ 
responding  to  Fr,p  =  0.  Hence,  it  Is  appropriate  to  make  the 
variables  dimensionless  in  the  following  way 

X  =  x'/T'jY  =  y'/T';N  =  T)'/T'jH  =  h'/T'jU  =  u'/(gT"  J1/2 

„  *1/2  ,  w  1/2  3/2 

V  *  v ' / ( gT '  )  ;F  =  <t>  +  i*  =  (♦'+  ir  )/g  T'  jP  =  p'/pgT'2 

C3-93 
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The  exact  boundary  conditions  of  two-dimensional  flow  are: 

2  ■'v 


U2  +  V2 


Pr, 


+  N  = 


T 


/  (Y  =  N(X) 


V  -  UN,X  =  0 


[3.10] 


[3.11] 


V  -  UH,X  *  0 


(y  =  h(x) )  [3.12] 


V  =  0,N  =  0,U  =  -  PrT  ( |  z|  -  »)  [3.13] 

W  =  U  -  iV  and  F  =  $  +  iY  being  analytical  functions  of 
Z  =  X  +  iY. 

(b)  Small  Perturbation  Expansion 

In  order  to  simplify  the  nonlinear  problem  we  seek  a 
solution  valid  for  small  Fr,p.  When  PrT  -*  0,  while  X,Y  =  0(l) 
the  flow  tends  to  rest  (Equation  3.13)  while  the  body  retains  its 
shape,  the  thickness  being  equal  to  unity. 

It  is  a  matter  of  simple  algebra  to  show  that  a  nontrivial 
small  perturbation  expansion  has  the  form 


U  =  PrT  Ui  +  PrT3U2  +  . . .  ^ 
V  =  FrTVi  +  FrT3V2  +  . . .  I 
N  =  PrT2Ni  +  PrT4N2  +  . . . 


[3.1*0 
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Tr.e  above  expansion  will  be  shown  to  be  regular  at  infinity 
and  consequently  there  is  no  need  to  consider  inner  and  outer 
expansions  separately. 

Expanding  U(X,Y)  and  V(X,Y)  in  the  vicinity  of  Y  =  N(X)  as 
given  by  Equations  [3.1^]  and  substituting  in  Equations  [3. 10]- 
[ 3 - 13 ] j  we  get  the  following  set  of  equations  after  separating 
terms  of  the  same  order 

(i)  Ui,  Vi,  Ni 


o 

II 

rl 

> 

(X  >  0,  Y  =  0) 

[3.15] 

Ni  =  |  (1  -  Ui2) 

(X  >  0,  Y  =  0) 

[3.16] 

Vi  -  UiH,x  =  0 

(Y  =  H(X) 

[3.17] 

Vi  ■  0;  Ui  -  -  1 

(|Z|  -  oo) 

[3.18] 

Hence  the  first  order  approximation 

the  free  surface  and  uniform  flow  at 

(ii)  Ua,  V2J  N2 

is  that  of  a  rigid 

infinity. 

wall  on 

V2  =  (UiNi),x 

(X  ,  -  0) 

[3.19] 

N2  =  -  UiU2 

(X  >  0) 

[3.20] 

v2  -  u2h,x  =  0 

(Y  =  H(X) ) 

[3. 2lJ 
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V2  =  0,  U2  =  0  ( fz|  -  « )  [3.22] 

In  the  second  order  approximation,  the  condition  on  the  un¬ 
perturbed  free  surface  (Equation  3-19)  is  equivalent  to  a  dis¬ 
tribution  of  sources  generated  by  the  first  order  flow,  with  no 
flow  at  infinity.  It  is  easy  to  ascertain  that 


[3.23] 


Since  Ni-*-  0  as  X  -*•  00  and  Ui  =  0  at  the  stagnation  point  X  =  0. 
Obviously  Wi  and  Wa  are  analytical  functions  of  Z. 

Higher  order  terms  satisfy  equations  similar  to  those  of 
second  order,  but  the  computations  become  tedious  as  the  order 
is  increased. 

c.  General  Solution 

The  solutions  of  the  different  order  approximations 
may  be  obtained  as  follows  (Figure  5):  The  region  AOBA  of  the 
Z  plane  is  mapped  on  the  £  half  plane  by 


Z  =  z(c)  [3.24] 

and  the  first  order  complex  potential  Fi  =  <3>i  +  i¥i  is  mapped 
on  the  same  £  plane  by 


Fi  *  const  X  £ 


[3.25] 
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At  second  (and  higher)  order  the  imaginary  part  of  i^2 

is  given  along  the  £  axis  (Figure  5b)  by  Equations  [3*19]  and 

[3.21] 


¥2  =  -  Uj.Ni  U  >  1) 

=  const  (i  <  l) 


[3.26] 


and  F2(C)  is  found  by  solving  the  related  Dirichlet  problem. 


(d )  Application  to  the  Rectangular  Body 

As  a  simple  example  we  consider  the  box-like  shape 
body  of  Figure  6.  The  AOBO  region  of  the  Z  plane  is  mapped  on 
the  £  half-plane  by 


Z  =  i  (Ca-  1)*  +  ~  in  [c  +  U2-  1)*]  [3.27] 


where  both  (Ca-  l)^  and  the  logarithm  have  real  determination  on 

i  >  1. 

With  Fi  »  -  Z/r  we  get 


Wi  «  Ui  -  iVi  -  %  (!~j  [3.28) 

In  particular  from  Equation  [3-16]  Ni  is  given  by 

Nx  -  \  (1  -  Uxa)  -  -5-  U  >  1)  [3.29] 
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Equations  [3»29]  and  [3.27]  describe  the  shape  of  the  free 
surface  in  a  parametric  form.  Ni(X)  is  represented  graphically 
in  Figure  7. 

The  next  order  term  F2  =  $2  +  i¥2  has,  according  to  Equa¬ 
tions  [3.21],  [3.28]  and  [3.20],  the  imaginary  part 

Y2  =  -  Uj.Ni  =  [|^J  U>  1,  n  -  0)  [3.30] 

¥2  =  0  ($<1,  n  -  0) 

and  W2  vanishes  at  infinity. 

F 2(C),  with  given  imaginary  part  on  the  real  axis  £  is 
determined  by  the  Cauchy  integral 

MO  -  -  f  MO  $  -  -  *  J  {5^372  $  [3‘311 

Y  1 

The  integration  in  Equation  [3.31]  may  be  carried  cut  analytically, 
the  result  being  for  C  ■  £,  £  >  1 

*<«>  ■  w  (iff +  (orf  l3-321 


The  velocity  Ua(0, O  is, 


accordingly 
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U2(UU)  *  ^  f  3  (— |xt"  '  j  +  fef  1 


(i+O*  \  ‘ 


I  £+1 


[3.33] 


and  finally  the  free  surface  elevation  is  given  by  Equation  [ 3 . 20 1 


No  =  -  Ui.Ua 


■'^-l  <  3  +  ■■  f‘g-r  .nU-(4?-  1)*] 

r^+i  >3  \  (r-  1)* 


[  3- 31*  ] 


N2  and  N/PrT2  =  Ni  +  Fr^N?  as  functions  of  X  are  represented 
graphically  in  Figure  7. 

( e )  Pressure  Distribution  and  Forces  Acting  on  the  Body 

The  dimensionless  pressure  has  the  following  expan¬ 
sion  resulting  from  the  Bernoulli  Equation, 


Fr 

P  -  -  Y  +  [l-(Ui*  +  V!*))-  FrT4(UiU8  *  VXV»  )+.... 


f3  351 

A  detailed  analysis  of  the  forces  acting  on  the  body  show  that 
the  drag  is  equal  to  zero,  as  it  should  be  in  an  ideal  fluid 
flow  with  no  waves.  The  dynamical  vertical  force  as  well  as 
the  moment  arc  different  from  zero  even  at  first  order.  The 
possibility  pf  computing  slnkage  an--,  trim  via  the  small  PrT  ex¬ 
pansion  will  te  explore*?  in  s  future  work. 
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(f )  Stability  of  the  Free  Surface 

Experiments  show  that  as  Pr^,  increases  a  breaking 
wave  appears  in  front  of  the  body  (Baba,  1969).  The  inspection 
of  Figure  7  reveals  that  as  Fr^,  increases  the  free  surface  be¬ 
comes  steeper.  Because  of  the  convexity  of  the  free-surface 
near  the  body,  the  centrifugal  effect  diminishes  the  pressure 
gradient  normal  to  the  surface.  When  the  pressure  gradient  be¬ 
comes  less  than  zero,  the  pressure  at  some  point  Inside  the  fluid 
is  smaller  than  the  atmospheric  pressure.  As  shown  by  Q. I. Taylor 
(1950),  such  8  condition  leads  to  the  instability  of  the  free- 
surface  and  very  often,  to  its  disruption  or  breaking.  Adopting 
the  vanishing  of  the  pressure  gradient  as  a  criterion  for  free- 
Burface  stability,  i.e.,  the  Taylor  stability  criterion,  we  are 
led  to  the  condition 


-  pgd  +  V,.*)"*  + 


u 


I*  +  vt» 


0 


Rewriting  Equation  (3.35)  in  dimensionless  variables  and 
with  r»  •  -  V  .  ,  (l+V»  1.  arrive  at  marginal  sta- 

|A  A 

blllty  for 

(U*  +  V^N^l+N^r3/2  -  -  (1+N,X#)’1/2  13.36) 

Expanding  Equation  (3.36)  on  the  free-surface  yields,  at  V 
order. 
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ST  *  '  +  aMW‘*xx> 

-1+|  Fr,p4Ni ,  *  +  PrT‘Nl>xNa,x  -  0  (3.371 

The  stability  criterion  as  expressed  by  Equation  [3.37]  has  been 
applied  to  the  flow  plast  a  box  shaped  body.  With  Mi,  N2,  U* 
and  U2  given  in  Equations  [3.28],  [3.29] ,  [3.33]  and  [ 3-3^1  the 
different  terms  of  [3*37]  have  been  computed  as  functions  of  §. 

In  Figure  7a  we  give  the  location  of  the  point  of  minimum  -dP/$n 
as  a  function  of  Fr^,.  The  point  of  minimum  -dP/dn  is  located  at 

X  «  0.3. 

(g)  Discussion  of  Results 

A  uniformly  valid  expression  for  the  velocity  and  free¬ 
st  face  profile  has  been  derived.  The  solution  of  first  order 
is  based  on  a  rigid  wall  approximation  while  in  the  second  order 
a  singularity  distribution  is  used  to  satisfy  the  free-surface 
condition.  The  solution  has  the  proper  behaviour  at  the  stagna¬ 
tion  point  S  (Figure  7)  since  both  W  and  dN/dX  vanish  there  In 
the  case  of  a  blunt  body.  The  behavior  at  Infinity  Is  also 
correct. 

Inspection  of  the  free-aurface  profile  as  a  function  of 
PrT(Pigure  7)  shows  that  as  PrT  Increases  the  free-aurface  be¬ 
comes  steeper.  This  Is  a  second  order  effect  and  reflects  the 
Influence  of  the  nonlinearity  of  the  free-surfaoe  condition. 
Although  at  Fr^  of  order  one  or  larger  it  is  doubtful  whether 
the  first  two  terms  represent  the  expansion  accurately,  the 
trend  Is  n  ;v?r*iheless  obvious. 
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The  pressure  gradient  normal  to  the  free-surface  decreases 
with  PrT  (Figure  8).  At  Fr^  =1.5,  Taylor  instability  of  the 
free  surface  occurs.  Of  course,  the  value  of  the  critical  FrT 
predicted  by  this  second  order  theory  is  probably  not  too  ac¬ 
curate,  but  the  calculation  serves  to  suggest  the  probable 
existence  of  a  critical  value  of  FrT,  beyond  which  wave  breaking 
occurs. 

In  analogy  with  progressive  free  surface  waves,  we  might 
even  expect  the  onset  of  Taylor  Instability  to  coincide  with 
the  non-existence  of  a  free  surface  wave  without  breaking. 

The  present  approach  permits  an  evaluation  of  the  influence 
of  the  bow  shape  on  the  inception  of  free  surface  breaking  as 
well  as  the  determination  of  the  ainkage  and  of  the  trim  of 
bodies  of  finite  length. 

A.  High  FJcoude  Number  Plow  (PrT  >  1):  The  Jet  Model 

(a)  General 

In  the  case  of  high  PrT  it  is  appropriate  to  relate 
the  variables  in  the  outer  sons  to  the  outer  length  U'*/g  and 
the  veloolty  U' 

x  ■  gx'/tJ'*;  y  *  gy‘/tJ**j  r\  -  gV/&‘*j  h  *  gh’/G’*; 

u  -  u*/G' i  v  -  v’/tt* j  a  -  g#ytt,s  (3.38) 
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The  exact  boundary  conditions  become  now 


U  +  V  .  1 

(y  «  n) 

[3.39] 

2  +  1  -  2 

u  -  v?},x  -  0 

(y  =  n) 

[3.40] 

u  -  vh,x  *  0 

(y  *  h) 

[3.41) 

u  -  1;  v  -*  0 

(x  -►» ;  y-*  -«) 

[3.42] 

with  w  ■  u  -  iv  an  analytical  function  of  z  «  x  t  iy. 

At  the  limit  T'g/tJ,a  -»■  (FrTa  -*  »)  the  body  collapses  into 
a  line  along  y  -  0  (Figure  3)  and  the  unperturbed  state  is  that 
of  uniform  flow.  The  first  order  equations  are  the  linearized 
equations  of  gravity  waves  of  the  type  [2.18]  -  [2.19]  (see  next 
paragraph). 

The  problem  of  two  dimensions,',  flow,  in  this  approximation, 
has  been  studied  extensively.  For  the  case  of  a  blunt  body  at 
the  free-surface  two  types  of  representations  have  been  suggested 
in  the  literature: 

(l)  Hie  replacement  of  the  body  by  a  source  (Wehausen 
and  Laltone,  I960)  or  by  a  constant  pressure  acting  on  the  free- 
surface  behind  the  bow  (Lunde,  1952).  It  is  easy  to  ascertain 
that  the  two  are  identical  If  the  source  is  located  on  y  ■  0. 

The  first  order  velocity  potential  for  a  source  of  strength 
Q  Is  (Vehausen  and  Laltone,  i960) 
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'  I  z 

f±(z)  =  4>i  +  ifi  =  ~  z  -  |  e‘iz  f  [3.^3] 

00 

The  free-surface  profile  corresponding  to  this  solution  has 
a  wavy  shape  far  behind  the  origin.  It  cannot,  therefore,  simu¬ 
late  a  semi-infinite  body  of  arbitrary  shape.  Near  the  origin, 
an  expansion  of  fi(z)  for  small  z  shows  that  the  free-surface 
is  continuous  there,  since  the  integral  in  [3*38]  behaves  like 
in  z  for  small  z.  The  complex  velocity  is  singular  near  the 
origin  like  in  z.  This  behavior  will  be  found  unsatisfactory 
for  matching  with  the  inner  solution  (paragraphs  c,  e). 

(ii)  The  replacement  of  the  body  by  a  pressure  dis¬ 
tribution  singular  at  the  leading  edge  like  xj_2  .  This  ap¬ 
proach  is  used  in  studies  of  planing  surfaces  (Sedov  1965, 

Maruo  1951,  Squire  1957)-  Approximate  solutions  have  been  found 
for  inclined  flat  plates  of  finite  length  by  a  Fourier  series 
expansion  of  the  pressure  distribution.  In  these  solutions  the 
velocity  w  is  singular  near  the  leading  edge  (z  =  0)  like  z  2, 
while  the  free-surface  is  continuous  there.  For  this  reason 
this  type  of  singularity,  although  stronger  than  that  of  (i), 

1 

is  still  too  weak  in  order  to  permit  matching  with  the  inner  ex¬ 
pansion  (paragraohs  c,  e).  An  interesting  feature  of  the  planing 
solution  is  the  fact  that  the  pressure  distribution  is  infe¬ 
rrable.  For  this  reason  the  leading  edge  correction  and  the 
inner  expansion  are  not  essential.  It  was  nevertheless  assumed 
that  a  jet  exists  at  the  leading  edge  and  Wagner  (1932)  has 
linked  the  jet  flow  and  the  pressure  singularity  in  a  way  simi¬ 
lar  to  the  matching  of  the  inner  and  outer  expansions. 
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Wu  (1967)  has  studied  the  flow  past  an  inclined  surface  in 
the  high  Pr  * egime  by  matched  asymptotic  expansion. 

x 

In  the  following  paragraphs  we  study  the  flow  past  a  blunt 
semi-infinite  body.  The  outer  expansions  corresponds  to  a  regime 
in  which  Fr  >  1  while  FrT  -*•  0.  Hence  we  are  in  the  range  of 

I  li 

displacement  ships ,  the  buoyancy  being  much  larger  than  the 
dynamic  lift,  dynamic  effects  being  important  only  in  the  bow 
region. 

(b)  The  Inner  and  Outer  Expansions 

The  inner  and  outer  expansions  of  Equations  [3-39]  - 
[3.42]  follow  closely  the  derivations  of  Section  II. 2,  the  body 
length  being  now  immaterial. 

The  outer  expansion  has  the  form  [2.16],  with  e  *  e*=]/ft»T2 
»  T'g/AJ'*  this  time.  Again  the  choice  of  the  first  order  ex¬ 
pansion  is  dictated  by  the  fact  that 

h(x)  =  c*H(x)  and  H(x)  =  0(l)  [3.44] 

The  first  order  terms  satisfy  equations  similar  to  [2.I7]  - 
[2.21]  which  may  be  written  in  a  complex  form  as 

Re  (wi  +  ifx)  =  0  (x  >  0,  y  =  0)  [3.45] 


=  T)i 


(x  >  0,  y  ■  0) 


[3.46] 
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^  =  H  (x  <  0,  y  =  0)  [3.^7] 

wi  =  0  (x  -*■  +  «)  [  3 . 48  ] 

fi  =  4>i  +  iYi  and  wi  =  ui  -  ivi  being  analytical  functions  of 
z  in  the  domain  y  <  0. 

The  inner  variables  are  those  of  [2.33]  with  Z  and  W  de¬ 
leted.  The  zero  order  inner  expansions  are  exactly  the  free- 
gravity  flow  Equations  [2.351  -  [2.38]. 

( c )  The  Zero  Order  Inner  Solution  for  the  Rectangular  Body 

Let  us  consider  again  the  simple  case  of  a  rectangular 
body  (Figure  9).  In  the  inner  limit  X  and  Y  are  fixed  with  re¬ 
spect  to  the  bow  and  we  assume  that  the  free-gravity  flow  there 
takes  the  form  of  a  jet  directed  upwards.  Gravity  effects  are 
taken  into  account  along  the  free-surface  upstream  by  the  outer 
expansion.  The  same  effect  on  the  jet  upwards  at  some  distance 
from  the  bow  is  ignored. 

The  solution  of  the  inner  problem  follows  the  classical 
methods  of  free  streamline  flow  studies  (Gurevich,  1965). 

The  complex  potential  Fq  plane  is  mapped  on  the  auxiliary 
£  =  £  +  ip.  half-plane  by 
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The  function  =  in  (l/WQ)  =  In  (l/Jw^j  )  +  10  has  given 
imaginary  and  real  parts  on  the  boundaries 


Ref 2  =  0 

o 

0  =  w/2 
0  =  3tt/2 

0  =  TT 


(A J ; t  >  o) 

(SJ;  -1  <  i  <  0) 
(SB;  -b2  <  i  <  -1) 
(BA;  i  <  -b2) 


)  [3-30] 


where 


b2  is  an  arbitrary  constant. 

The  mapping  of  on  (  is  a  solution  of  a  classical  mixed 

problem  (Signorinl  problem)  which  is  reduced  a  Dirichlet 

l 

problem  for  the  function  Qq/Fqp .  The  result  of  the  integration 
of  Cauchy's  integrals  is 


Q.  =  in 
o 


1 

1  ^ 

r 

2 

£*  -  i 

C  +  ibi 

+  i 

U*  -  ibj 

+  iv 


[3.51] 


With  a  new  auxiliary  plane  o  related  to  £  through 


0)  =  £ 


[3.52] 


we  get  from  [3-51]  and  [3.52] 


W 


0 


0+1 

o-i 


p-lb\^ 
o+ib  | 


[3.53] 
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Eigensolutions  of  the  type  ~  j_£n+l/2  are  ru]_e(j  out  since  they 
yield  infinite  or  zero  velocity  far  downstream. 


The  mapping  of  the  Z  plane  on  the  co  plane  results  from  the 
basic  relationship 

dF 

dZ  -  -jp  [3.54] 

0 


Substituting  [ 3 - ^9] >  [3*52]  and  [3.53]  into  Equation  [3-54] 
and  taking  in  consideration  that  Z  =  -  i  for  a)  =  -  ib  (Figure  8) 
we  get 


Z  +  i 


2A  j  afH  cD-i  co+lb  | 2 

ir  I  a)  CD+i  co-ib  j 


-ib 


[3.55] 


The  integral  of  [3-55]  can  be  carried  out  in  a  closed  form  with 
the  result 


Z  +  i  =  f  /  |  (o>2+  b2P  -  %-  in  ^-+  ~  -  i (2-b ) (o)2+b2 


-ib 


v. 

+  (gb.!)  tn  +  !  tn  t>  ±M±  b2)2  1) 

'  -ib  ico  1 


[3.56] 


In  all  the  above  formulae  the  square  roots  and  the  loga¬ 
rithm  have  real  determination  on  <jo  =  real. 


The  unknown  constants  A  and  b  have  to  be  determined  from 
matching  with  the  outer  expansion.  For  this  purpose  let  us 
seek  the  behavior  of  Z  and  Wq  far  from  the  bow,  i.e.  for 
| od |  »  b.  From  Equation  [3*56]  we  get  for  large  o 
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Z  +  I  .  ^  f 

7T  2 


-  i(2-b)  cu  +  (2b-l-~-  )  in  03  +  i  \  (2b-l-^-  ) 


+  + 

1  1  •  •  • 

03 


[3.57] 


Similarly,  from  Equation  [ 3 • 53 ]  we  obtain  for  W 


W  =  -  1  - 


i(b3- 


+...  [3.58] 


Two  cases  of  interest  are  to  be  discussed  separately: 

(i)  b  ^  2.  In  this  case  in  a  first  approximation 


W  =  -  !  -  .ug=£tei  1  + 


[3.59] 


Nq  -  -  1  +  (2b-l-^  )A-2(2-b)'^|  X^+ - (X  >  0)  [3.1 


No  *  -  1 


(X  <  0)  [3.61] 


Hence,  the  velocity  perturbation  behaves  like  Z  ,  while 

N  ~  -  X*. 
o 

(li)  b  *  2.  For  this  distinguished  value 


[3.62] 
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4a3/2 

Nc  =  -  1  +  A  +  (X  >  0)  [3.63] 

N  =  -  1  (X  <  0)  [3.64] 

0 

In  this  case  for  large  X  the  velocity  perturbation  decays  like 

_o/2  -1/2 

Z  J/  .  while  N  tends  to  the  constant  value  -1  +  A  like  X  '  . 

*  o 

There  is  no  major  difficulty  in  determining  the  zero  inner 
solution  for  bodies  of  other  shapes  than  the  rectangular,  pro¬ 
vided  that  0  is  given  as  a  function  of  C  (for  instance,  a  poly¬ 
gonal  body).  In  the  case  of  an  arbitrary  body  with  given  0  as 
a  function  of  x  (or  y)  the  problem  becomes  extremely  difficult 
and  leads  to  an  integral  equation  for  0(C)  (Wu,  1967). 

(d)  The  First  Order  Outer  Solution 

The  outer  problem  reduces  to  the  determination  of 
fi(z)  subject  to  conditions  [3.451*  13.46]  and  [3.46],  while 
for  the  particular  case  of  a  rectangular  body  Equation  [3.471 
gives 

¥1=1  (x  <  0,  y  -  0)  [3.65) 

The  problem  is  made  unique  if  the  singular  behavior  of 
fi  (or  wi)  near  the  bow  (z  -  0)  is  prescribed.  The  inner  so¬ 
lution  sh  )ws  that  there  exist  two  possibilities  for  WQ:  Equa¬ 
tion  (3.591  or  Equation  [3.62].  A  detailed  study  shows  that 
matching  is  possible  only  in  the  second  case.  The  reason  is  the 
following:  if  b  /  2  the  matching  requires  that  wi  ~  z“'  near 
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the  origin  and  A  or  (2-b)  has  to  be  of  the  order  In  the 

first  case  qi  is  continuous  at  x  =  0  and  has  the  value  qi  =  -  1 
there;  this  requires  a  solution  with  tjx  dropping  from  qx=0(x  =  «) 
to  tjx  =  -  1  (x  =  0).  Such  a  solution  is  not  possible  for  the 
assumed  type  of  singularity  of  wi .  If  (2-b)  is  0(€*^)  and  A  =  0(l), 
Tjx  has  a  jump  at  x  *  0  from  (-1  +  A)  to  (-1).  Again  the  assumed 
type  of  singularity  of  wi  does  not  allow  for  a  discontinuous  tjx 
(see  paragraph  a  of  this  section). 

Consequently,  we  adopt  the  value  b  -  2,  and  the  inner  term 

contains  A  as  the  only  unknown.  Moreover,  wi  behaves  near  z  =  0 

-V2  -4 

like  z  '  while  tjx  is  singular  like  x  for  x  >  0. 

An  exact  solution  of  fi(z)  is  still  difficult.  The  usual  way 
to  find  it  (Sedov,  1965 )  is  to  consider  the  function  wi  +  ifx 
(suggested  apparently  for  the  first  time  by  Keldlsh)  and  to  con¬ 
tinue  it  analytically  over  x  >  0  in  the  entire  z  plane  cut  along 
x  <  0.  With 

wx(z)  +  ifx  (z )  -  k(z)  (3.66) 


the  unknown  function  k(z)  has  to  be  imaginary  for  z  ■  x  >  0.  Its 
real  part  is  in  fact  the  linearized  pressure.  The  solution  of 
Equation  [3.66]  with  the  radiation  condition  [3.48]  Is 


fi(s) 


e1Xk(*)dX 


13.67] 
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The  function  r(-£,  -iz)  is  analytical  in  the  whole  z  plane  cut  by 
x  *  0,  y  >  0.  It  has  the  following  asymptotic  series  (Qradshtein 
et  al,  1965): 

(i)  For  small  z: 

r  00  1 

j  _  /ii  1  I  * —  /  ,  \n  #  .  _  \n  I 


r(-i,  -iz)  -  r(-i)  +  elir//J|  z~*  2  -Y"  l-.1!  (~lz) 

L~.  n.'  (n  -  i) 
n-1 


[3.72] 


Hence  with  r(-£)  *  -  2 r*,  fi  has  the  expansion 


fx(z)  -  -  12tt*  a  e*1^  +  2ia  +  0(z^  C 3 . 73 ] 


Yx(x,0)  -  -(2*)*  a  +  2ax“*  +  0(x*)  (x  >  0) 

I 

¥x (x,0)  -  -(2r)*  a  +  0(  x  *)  (x  <  0) 


[3.74) 


(ii)  For  large  z: 


r(-i,  -iz)  -  |,|'3/8  .-13(,re  .l*fi  ♦  oj-jij-j 


fi  ■  -a  e 


■13  zrg  z/2  H  +  0|_i_j 


13-75) 


[3.76) 
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Agaln,  Yi  *  Imfi  follows  the  expressions 


(x  >  0) 
(x  <  0) 


[3.77] 


Unfortunately  Yi  is  not  constant  along  x  <  0,  as  required 

by  Equation  [3.65].  But  the  approximate  solution  has  the  proper 

-3/2 

behavior  near  the  bow,  where  wi  ~  z  '  and  *  Yj  is  like 

for  x  >  0,  and  also  at  x-*  •  •  where  tjx  —  -  a  with  no  waves 
left  behind  the  body. 

Now,  It  is  a  matter  of  convention  how  to  pick  the  value  of 
a  in  order  to  satisfy  approximately  Equation  [3.65].  If  we  try 
to  satisfy  [3. 56]  near  the  bow  a  may  be  obtained  from  the  crn- 
d i tlon 

fi  (0,-0)  -  1  (3.78) 

which  gives 

a  -  -  1/(2*)*  (3.791 


Although  we  have  no  exact  solution  for  the  outer  problem, the 
approximate  expression  f 3.71]  reflects  the  main  features  of  the 
solution. 

(f )  The  Matching  of  the  Inner  and  the  Outer  Solutions 

The  matching  Is  generally  achieved  by  an  intermediate 
expansion  (Cole,  1968).  In  the  present  case  it  can  be  done  by 
the  simple  principle  (Van  Dyke,  196*):  The  outer  limit  of  the 
Inner  solution  equals  the  inner  limit  of  the  outer  solution. 
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Subs  titu  ting  z  =  e*Z  in  the  outer  solution  [3*71]  and  seeking 
the  limit  Wi(Z)  for  €#  -*  0  and  Z  *  0(l)  we  obtain  from  Equa¬ 
tions  [2.16]  and  [3.73]  for  the  inner  limit  of  the  outer  solution 


3 75*375  -  1  +  °"*1/2) 


[3.80] 


The  outer  solution  matches  with  the  outer  limit  of  the  inner 
solution  [3*62]  only  if 


„  *-1/3 


A  ~  e 


[3.81] 


The  estimate  of  [3.8l]  is  the  main  result  of  our  analysis. 
In  particular,  for  the  value  of  a  of  [3*79h 


.  ,-1/3 


^  *  2  6 


.-1/2 


13.82] 


The  matching  of  tj  and  NQ  is  also  ensured  at  order  c  '  with 
A  given  by  [3.82]:  Prom  Equations  [2.16)  and  [3.7*0  we  find  for 
the  inner  limit  of  the  outer  solution 


H  ■  - 


+  0(1) 


13.83) 


€**  X* 


while  the  outer  limit  of  the  inner  solution  [3.83]  has  the  form 


€*?  XJ 


(3.84) 
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If  we  assume  that  the  bow  drag  for  a  body  of  finite  length  L' 
has  the  same  expression  we  have  for  the  bow  drag  in  its  con¬ 
ventional  form 

D'/pTJ,2L'  =  (T'/L»  )(U'2/ST'  )T  =  (T'/L'  )T(Hi'2/gL'  r  [3-86] 

(h)  Discussion  of  Results 

In  the  present  section  the  free-surface  flow  past  a 
blunt  body  with  high  PrT  number  (but  low  Fr^)  has  been  studied. 
The  problem  is  different  from  that  considered  in  planing  studies 
since  the  position  of  the  body  is  fixed  and  its  bottom  is  hori¬ 
zontal. 
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The  main  results  of  the  analysis  are  the  following: 

(i)  The  proper  type  of  pressure  singularity  at  the 
bow  in  the  outer  solution  is  of  the  order  |x|~^2.  This  pressure 
is  not  integrable  so  that  lift  may  be  evaluated  only  via  the 
inner  expansion.  Obviously,  the  inner  solution  shows  that  the 
dynamic  pressure  is  a  maximum  pU,2/2  at  the  stagnation  point. 

(ii)  A  Jet  is  assumed  to  appear  at  high  Pr,p  numbers. 

The  energy  of  the  jet  is  probably  entirely  dissipated.  This  is 
causing  a  drag  additional  to  the  wave  resistance. 

(iii)  The  jet  thickness  and  the  bow  drag  grow  slowly 

with  FrT,  like  . 

The  present  analysis  may  be  refined  in  different  directions: 
By  improving  the  outer  solution,  by  considering  bodies  of  finite 
length,  by  studying  different  bow  shapes  and  by  extending  the 
results  to  real  flat  ships. 
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IV.  CONCLUSIONS 

Special  approximations  are  needed  in  order  to  analyze  the 
free  surface  flow  in  the  vicinity  of  the  bow  of  blunt  ships.  In 
the  case  of  thin  ships  (T'/B'  sufficiently  small),  the  inner  bow 
flow  reduces  to  a  two-dimensional  flow  in  a  vertical  plane  normal 
to  the  bow.  Furthermore,  it  is  appropriate  to  consider  a  blunt 
two-dimensional  body  of  semi-infinite  length  and  this  is  done 
herein. 

The  situation  at  the  confluence  of  a  blunt  bow  and  the 
free  surface  is  clarified  first.  It  is  shown  that  there  are  two 
possible  angles  between  a  free  surface  and  a  rigid  wall  at  the 
stagnation  point:  (i)  if  the  wall  is  inclined  with  respect  to 
the  horizontal  at  an  angle  larger  than  120°,  the  free  surface 
intersects  the  wall  at  120°;  (ii)  if  the  wall  is  inclined  at 
less  than  120°,  the  free  surface  is  horizontal.  The  latter  case 
is  usual  for  ships. 

For  small  Froude  number  based  on  draft,  Fr^.,  the  flow  can 
be  analyzed  by  means  of  an  expansion  in  FrT,  according  to  which 
the  first  approximation  corresponds  to  replacing  the  free  sur¬ 
face  by  a  rigid  wall.  The  flow  past  a  rectangular  body  is 
analyzed  to  the  second  order.  The  solution  to  this  second  ap¬ 
proximation  shows  that  the  free  surface  becomes  steeper  in 
front  of  the  bow.  Application  of  Taylors  instability  criterion 
leads  to  the  conclusion  that  the  stability  of  the  free  surface 
decreases  with  increasing  Fr^,  presumably,  culminating  in  breaking. 
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In  the  second  approximation,  the  free  surface  becomes  unstable 
at  PrT  of  about  1.5  and  at  a  position  30  percent  of  the  draft 
ahead  of  the  bow.  The  small  Pr^,  theory  does  not  allow  the  cal¬ 
culation  of  bow  drag,  which  only  ensues  after  breaking,  but  it 
does  permit  the  estimation  of  sinkage  and  trim. 

For  large  FrT  (but  small  Fr^)  the  flow  past  a  rectangular 
bow  has  been  analyzed.  The  problem  is  different  from  that  con¬ 
sidered  in  planing  studies,  since  the  bow  is  vertical,  while 
the  bottom  is  horizontal.  The  problem  is  solved  by  matching 
appropriate  inner  and  outer  solutions.  The  inner  solution  cor¬ 
responds  to  a  free  surface  without  gravity  while  the  outer  flow 
corresponds  to  the  usual  linearized  free  surface  flow  with 
gravity.  The  main  results  of  the  analysis  are:  (i)  the  proper 
type  of  pressure  singularity  at  the  bow  in  the  outer  solution 
if  of  order  |x|^2;  (ii)  a  spray  jet  appears  at  the  bow,  whose 
energy  is  probably  entirely  dissipated.  This  jet  causes  a  bow 
drag  additional  to  the  usual  wave  resistance;  (iii)  the  jet 
thickness  and  the  bow  drag  grows  slowly  with  Fr,p,  like  Fr,^1/^. 
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FIGURE  2  -  BEAM  DRAFT  VERSUS  (FROUDE  NUMBER )0IUFT 


FIGURE  7  -  THE  FREE -SURF ACE  SHAPE  IN  FRONT  OF  A  RECTANGULAR  BODY 
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FIGURE  8  -  THE  RELATIONSHIP  IE  TWEEN  THE  MINIMUM  PRESSURE  GRADIENT 
NORMAL  TO  THE  FREE -SURF  ACE  AND  THE  FROUOE  NUMKR 


HYDRONAUTICS,  INCORPORATED 


FIGURE  9  -  FREE- GRAVITY  JIT  FLOW  PAST  A  RECTANGULAR  |<H)Y 
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- Special  approximations  are  needed  in  order  to  analyze  the 

free  surface  flow  in  the  vicinity  of  the  bow  of  blunt  ships.  In  the 
case  of  thin  ships  the  inner  bow  flow  reduces  to  a  two-dimensional 
flow  in  a  vertical  plane  normal  to  the  bow.  The  situation  at  the  con¬ 
fluence  of  a  blunt  bow  and  the  free  surface  is  clarified  first.  It 
is  shown  that  there  are  two  possible  angles  between  a  free  surface 
and  a  rigid  wall  at  the  stagnation  point:  (i)  if  the  wall  is  inclined 
with  respect  to  the  horizontal  at  an  angle  larger  than  120°,  the  free 
surface  intersects  the  wall  at  120°;  (fl)  if  the  wall  is  inclined  at 
less  than  120°,  the  free  surface  is  horizontal.  The  latter  case  is 
usual  for  ships . ;  -  s 

For  small  Froude  number  based  on.  draft,  FrT,  the  flow  can  be 
analyzed  by  means  of  an  expansion  in  Frij,  according  to  which  the  firs 
approximation  corresponds  to  replacing  the  free  surface  by  a  rigid 
wall.  The  flow  past  a  rectangular  body  is  analyzed  to  the  second  order; 
Application  of  Taylors  Instability  criterion  leads  to  the  conclusion 
that  the  stability  of  the  free  surface  decreases  with  increasing  FrT, 
presumably,  culminating  in  breaking.  The  small  FrT  theory  does  not 
allow  the  calculation  of  bow  drag,  which  only  ensues  after  breaking, 
but  It  does  permit  the  estimation  of  sinkage  and  trim. 
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